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SUMMARY 


Multilevel self-adaptive Newton-Raphson type strategies are developed to 
improve the solution efficiency of nonlinear finite element simulations of 
statically loaded structures. The overall strategy involves three basic levels. 
The first level involves preliminary solution "tunneling" via primative opera- 
tors. Secondly, the solution is constantly monitored via so-called quality/ 
convergence/nonlinearity tests. Lastly, the third level involves self-adaptive 
algorithmic update procedures aimed at improving the convergence characteris- 
tics of the Newton-Raphson strategy. Numerical experiments are included to 
illustrate the results of the procedure. 


INTRODUCTION 

Finite element (FE) or difference simulations of continuum problems gener- 
ally lead to nonlinear modelling equations [1,2]. Generally, such simulations 
must be solved by various techniques which are inherently iterative in nature. 
For instance, such methodologies as direct numerical integration, Newton-Raph- 
son (NR), and modified Newton-Raphson (MNR), as well as the incremental ver- 
sions of such procedures (INR, IMNR) have all been employed [2]. Since the 
types of nonlinearity exhibited by continuum problems are both diverse and com- 
plex, the question of the best choice of an appropriate solution algorithm 
inevitably arises. Note, while many alternatives are available, generally the 
various solution procedures may have special advantages for certain classes of 
problems but may exhibit poor convergence for other situations. 

In this context, the ideal general purpose (GP) nonlinear FE code should 
have numerous algorithmic options augmented with a degree of artificial intel- 
ligence. Namely, the problem solving capability should involve a heuristi- 
cally guided trial and error search in the space of possible solution via an 
automatically structured algorithm. Unfortunately, because of the inherent 
difficulties associated with code architecture and kinematic, kinetic, con- 
stitutive and boundary condition formulations, generally only one algorithmic 
option is usually available in GP codes. In this context, because of its wide 
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applicability, most GP nonlinear FE codes employ some variant of either the 
straight or modified INR algorithmic procedures. 

Note, while nonlinear codes present the user with far reaching capabilities, 
without a priori physical insight, expensive parametric studies are oftentimes 
necessary to insure adequate solution convergence. For instance, unless the 
proper load increment is employed, either poor convergence or out of balance 
loads are generally encountered. Incorporating heuristic programming could 
eliminate some of the expensive and time consuming parametric studies that are 
now required to determine the proper incrementation necessary for reasonable 
convergence. 

In view of the shortcomings of the current generation of solution al- 
gorithms, this paper will consider the development of self-adaptive NR stra- 
tegies for the solution of nonlinear FE or difference simulations of stati- 
cally loaded structures. The main thrust will be to consider strategies which 
for the most part are compatible with currently available GP codes. The over- 
all development will be considered in three main levels. The first will in- 
volve the use of INR operators to "tunnel" into the solution space in the usual 
manner. The second level will involve the constant monitoring of the different 
stages of solution via various qual ity/convergence/nonl inearity tests. Finally, 
the last level is an outgrowth of the findings of the second; namely, if one or 
more of the qual ity/convergence/nonl inearity tests are violated, various sce- 
narios are then triggered to modify the INR strategy. 

Based on the foregoing, the paper will outline in detail the multilevel 
static solution strategy, the development of the qual ity/convergence/nonl inea- 
rity tests as well as overview the various self-adaptive iterative update pro- 
cedures. The analytical considerations will be complemented by several nume- 
rical experiments which outline the various aspects of the quality/convergence/ 
nonlinearity tests and which demonstrate the self-adaptive strategy. 


MULTILEVEL SOLUTION STRATEGY: 
OVERVIEW 


As noted earlier, unless the proper load incrementation is employed, 
either poor convergence or out of balance loads are generally encountered. 
Such anomalous behavior is generic to all nonlinear codes employing non-self- 
adaptive INR algorithms. In this context, the main thrust of this work is to 
establish a three level iterative solution strategy involving: 

i) Level 1; Preliminary solution development via the primative but 
computationally efficient IMNR algorithm; 

ii) Level 2; Solution monitoring via qual ity/convergence/nonl inearity 
tests and; 

iii) Level 3; Self-adaptive update procedures to modify the primative 
operator. 
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Note here computational efficiency is meant to be a measure of the amount of 
time spent during a cycle of iteration not the overall process. 

The main purpose of the first level of the overall strategy is essentially 
twofold. The first is to generate the most efficient solution if the requisite 
qual ity/convergence/nonl inearity criteria are satisfied. If not, the infor- 
mation generated by the IMNR "tunneling" of the solution space can, through the 
second level tests, trigger the proper third level action. 

In terms of the foregoing, it follows that the second level is essentially 
threefold in nature. The quality check involves monitoring: the rate of con- 

vergence; monotonicity; positive, negative and semi-definiteness; etc. The con- 
vergence tests check for outright solution failure, and lastly, the nonlin- 
earity tests ascertain the "degree" of nonlinearity excited. 

In the third level, the foregoing information is used to trigger various 
self-adaptive modifications of the IMNR iterative strategy. Namely: 

i) Global stiffness reformation; 

ii) Preferential local reformation and; 

iii) Load increment adjustment. 

Such algorithmic adjustments form the heart of the third level of the overall 
strategy . 


INR FAMILY OF STRATEGIES 


The overall family of INR strategies can essentially be established by 
introducing increasingly severe restrictions to the straight methodology. 
Specifically, starting with the virtual work theorem depicted by [1] 
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the typical FE shape function formulation yields the following nonlinear large 
deformation field equations [1] 
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such that g, S, Y and F are, respectively, the strain tensor in vector form, 
the second PiolaKirchhoff pseudo stress tensor, the nodal displacement vector, 
and lastly, the nodal force vector. 

To solve (2), the Taylor expansion theorem can be used to establish the 
following tangent stiffness formulation, namely 

[Kt(Yj_j)] AyJ = AF. (5) 

where AY^ denotes the i^*^ nodal displacement iterate associated with the 
load increment. The nodal displacement, tangent stiffness and load imbalance 
are defined by 
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such that I|(, [a(Y.j_J], [D^(Y.j_j)] and F^ respectively denote the number of 
iterations required of Toad step, the initial stress matrix, the tangent 

material stiffness and the total nodal load after increments. For the strai- 
ght INR approach, [K^] is continuously reformed and inverted. This is obvious- 
ly quite expensive. In this context, the following versions of the INR algori- 
thm can be established for a specific load increment solution cycle, namely; 


i) 

ii) 

iii) 

iv) 

V) 

vi) 
vi i) 


Straight INR with constant reformation of tangent stiffness matrix 
during iteration; 

Intermittent global reformation during iteration; 

Preferential local reformation during iteration; 

BFGS type [3] reformation during iteration; 

Classical modified INR procedure wherein stiffness is reformed only 
at beginning of load step; 

No reformation, just iteration; 

Reformation with no iteration, etc. 


As will be seen later, various versions of the foregoing INR family of al- 
gorithm are incorporated in the self-adaptive strategy. This will obviously 
lead to a hierarchy with varying degrees of computation power/efficiency. 
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QUALITY/ CONVERGENCE/ NONL I NEARI TY TESTS 


The qua! ity /convergence /non! inean'ty tests are the core of the multi- 
level strategy. Such tests are themselves organized into three main cate- 
gories, namely: 

i) Classical norm type convergence tests; 

ii) Quality of convergence tests and; 

iii) Degree of nonlinearity tests. 

The first group of tests are essentially of the normed type pass or fail vari- 
ety as typified by: 

a) The out of balance norm test; 

l|4F!-i|/||AF'.J|<tol (9) 

b) The global displacement norm test; 

II yJII /II Y^_J| < tol (10) 

The main intent of such tests is essentially to monitor the success or failure 
of the iterative process. Note, while such tests are efficient and well adap- 
ted to this purpose, they cannot be used effectively to forecast potential 
difficulties until outright failure occurs. 

In this context, what is required are so-called quality checks which en- 
able a constant monitoring of the solution so as to determine whether the direc 
tion of convergence is proper. This is the purpose of the second stage of 
checking. Namely, the quality checks test whether the iterative process pos- 
sesses the requisite: rate; monotonicity; positive, negative and semi definite 

ness; etc. Once determined, such information is used to trigger the various 
modifications of the primative first level IMNR strategy. 

Since the paper is mainly concerned with static loading problems, various 
statements concerning the quality of solution convergence can be made at the 
outset. For instance, since most static loading is applied in a monotone fash- 
ion, it is expected that unless there is overshoot, successive iterated 
solutions should behave as a monotone, positive, negative or semi definite se- 
quence. Behavior to the contrary obviously represents either overshoot or po- 
tential divergence. 

Since it is difficult to ascertain the monotonicity and definiteness 
from either of the normed or vectorial versions of the nodal displacements and 
forces, alternative field measures must be employed. In this direction, the 
local (element) and global strain energy stored can serve in such a capacity. 
This follows from the fact that for monotone loading situations, successive 
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iterations lead to a monotone positive definite sequence of energy iterates 
for softening structure. In the case of hardening situations, successive ite- 
rates may be nonmonotone for at least the first two iterates. Thereafter, the 
energy iterates tend to be monotone and negative definite. This process is 
clearly seen by the normed analogy of the iterative process depicted in Figures 
1 and 2. 

The incremental iterate energy stored during a given iteration step is 
essentially the shaded area illustrated in Figure 1. Realizing that the or- 
dinate values of the true solution curve are given by (1), it follows that the 
incremental energy stored during the iteration step of the load incre- 
ment can be approximated by the following inner product, that is 
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Assuming that a total of K iteration steps are associated with the a load- 
step, then the following expression can be developed for the energy stored, 
namely: 
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Note (12) is essentially a trapezoidal type integration approximation for 
the area under the hyper-curve defining the solution of the loadstep. Now, 
summing (12) over the entire set of L loadsteps associated with a given problem 
yields the requisite overall strain energy stored, namely : 
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To obtain the strain energies for say the e^^ element, (13) must be inter- 
preted from a local point of view. Namely, the requisite partitions of F|^® 
and AY|^® must be employed in a partitioned version of (13), that is 


Cal-lj, 


(14) 


where here AY and F|^ , are respectively the local and element nodal displace- 
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merits and forces. Note, due to the form of (14), any form of tangent stiffness 
type of constitutive law can be accomodated. 

In terms of the iteration process associated with softening media, it fol- 
lows that for convergent situations 

> ...> ...> 0 (15) 

that is, successive iterates are monotone and positive definite. Hence, for 
softening media, a check of successive iterates for monotone decreasing posi- 
tive definiteness will establish a measure of the quality of convergence. For 
the hardening case, the E^ sequence associated with the convergent solution 
process takes the form: 

E^ < E^ <...< Ef < eJ., ...<0<E^ (16) 

2 3 K K"t| X 

0 

As can be seen, for k>l, the sequence E|^ is monotone increasing but negative 
definite. In this context, similar to the softening case, a check of the mono- 
tone increasing negative definiteness of successive iterates is used as one of 
the measures of the quality of convergence. 

A last but very important way of predicting potential solution difficul- 
ties can be achieved by monitoring the degree of nonlinearity excited as the 
deformation process continues. This can be achieved by selectively checking 
the changes of curvature of the global and local strain energy space. Such 
behavior can be ascertained using difference operators to evaluate either the 
slope, rate of change of slope, or more elaborately, the radius of curvature of 
the energy space as either a function of the loading parameter or the nodal 
displacements. An alternative approach would be to locally spline fit the ene- 
rgy-loading parameter space. In this way, the current curvature/slope can be 
obtained either on a local element or global basis. Such information can be 
used to initiate changes in load step size as well as to control local and 
global stiffness reformation. 

The importance of such tests follows from the fact that although FE simu- 
lations of structures composed of general media undergoing large deflections 
are inherently nonlinear, the degree of nonlinearity excited varies from point 
to point as well as from load increment to load increment. As it is possible 
that large portions of the structure may exhibit basically linear behavior, 
many general purpose codes allow the user to partition the overall structure 
into its linear and nonlinear groups. Although this certainly adds to the ef- 
ficiency of the code, generally such information is not known a priori unless 
extensive parametric studies have already been performed. In this context, the 
nonlinearity check will enable the automatic partitioning of the structure by 
allowing for preferential reformation of the tangent stiffness depending on the 
amount of local nonlinearity excited. 
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ADAPTIVE STRATEGY 


In the context of the Inherent features of the INR family of algorithms, 
the adaptive strategy incorporates the following procedural options namely: 

i) Tangent stiffness reformation and; 

ii) Load increment adjustment. 

Each of these options in turn involves several different levels. For instance, 
stiffness reformation can be considered in several stages, that is: 

i) Global reformation; 

ii) Preferential local reformation or; 

iii) 6FGS [3] reformation. 

The adaptive incremental load can also be achieved in several ways namely: 
i) Increment expansion; 
ii) Increment contraction or; 

iii) Corrective incrementation. 

As noted earlier, the initiation of either option is dependent on three basic 
criteria, that is: 

a) Quality of convergence; 

b) Outright failure to converge or; 

c) The degree of nonlinearity excited. 

While the reformation option is triggered by the second level tests, the 
specific adaption triggered is primarily dependent on the degree of nonlinear- 
ity excited. Hence, for mildly nonlinear (elastic) situations, the BFGS refor- 
mation process is employed. In the case where significant local or global non- 
linearity is excited, then either global or preferential reformation is ini- 
ti ated. 

As can be seen from the proceeding categories, various types of load in- 
crementation are possible. The overall strategy is a combination of such 
options. Specifically, when significant solution degradation is monitored by 
the level two tests, then corrective incrementation is initiated. Namely, ne- 
gative load incrementation is employed to enable the retracing of a portion of 
load history wherein a lower order algorithmic strategy yielded poorly conver- 
ged results. 

To strike a balance between solution convergence and economy, the overall 
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adaptive strategy is centered about a primative version of the INR algorithm 
namely the IMNR. Depending on the results of the qual ity/convergence/nonl inea- 
rity tests, the level of the IMNR is either upgraded or lowered by modifying 
the pattern of stiffness reformation and incrementation. Note, since the main 
incentive is to achieve a successful solution at least cost, the hierarchy is 
ordered to first implement increment adjustment and then reformation. As a 
further move to achieve economy, the global reformation process typically em- 
ployed at the start of an IMNR increment can be established preferentially de- 
pending on the curvature tolerance associated with the global and local non- 
linearity checks. 

Since space is limited, a full description of the various detailed hier- 
achies must be left to future publications. In this context, for the present 
purposes. Figure 3 gives a good overview of all the possible flows of control 
associated with the three-level strategy. As can be seen from this figure, 
contingent upon the various "flags" generated in the level two tests, the con- 
dition code check routine will initiate the actual modification of the INR 
strategy along the lines outlined in the proceeding discussion. 


DISCUSSION 

Interestingly, while such factors as geometry, material properties, 
boundary conditions, etc., all have some effect on the choice of load incre- 
ment size, once an excessive value has been chosen, typically similar types 
of solution degradation are encountered when only the primative non-self-adap- 
tive algorithm is used. Specifically, three basic types of solution pathology 
tend to occur. These can be categorized by: 

i) Immediate and strong nonmonotonicity; 

ii) Moderate but progressively increasing nonmonotonicity and non- 
positive definiteness and; 

iii) Mild monotonicity with either very gradual increases or decreases 
in solution oscillation. 

Note, such behavior can be excited either in the first or successive load steps. 
Figures 4 and 5 give examples of such behavior. While the results illustrated 
pertain to a rubber sheet, similar results were obtained for elastic/plastic 
media as well as for different geometries and boundary conditions. 

The solution failure depicted in Figure 5 is typical of those that usually 
arise. Specifically, as can be seen, for the given load increment excellent 
convergence is obtained in the first step. In the second, a mild form of non- 
monotonicity and nonpositive definiteness is encountered. Finally, in the 
third step, strong and progressively increasing nonpositive definiteness is 
encountered. Solution failure is finally initiated by out of balance loads. 
This scenario is typical of excessive load incrementation. Note, as can be 
seen from these results; the onset of such behavior is signalled by the ini- 
tiation of nonmonotonicity or incorrect definiteness. By studying the behavior 
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of local element energies, additional insights are obtained. For the problems 
illustrated in Figures 4 and 5, the solution degradation is initially loca- 
lized but gradually spreads to the entire structure as the iteration process 
continues. Employing the self-adaptive strategy to the foregoing problems 
caused the second level monotonicity tests to trigger automatic increment ad- 
justment and preferential stiffness reformation. This led to the generation 
of the correct solution. The overall strategy was tested on several nonlinear 
problems which exhibited pathological behavior for given load increment choices. 
The types of problems considered combined varying degrees of kinematic, kinetic 
and material nonlinearity. In each case barring possible bifurcations, the 
level two tests were able to automatically initiate the requisite corrective 
self adaptions to enable successful solutions. The .main problem encountered 
with the concept of self-adaptive strategies arises from the fact that some 
engineering insight must be practiced in order to cut down overall running 
times. Otherwise, excessive execution times are encountered as the adaptive 
strategy shifts gears to adjust for improper incrementation. 
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Figure 1.- Iteration process for hardening problem. 



Figure 2.- Iteration process for softening problem. 
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GLOBAL ENERGY INCREMENT (iNCH-LBS.) 



Figure 3.- Overall flow of control of three-level self-adaptive strategy 



2 3 4 

Figure 4.- Global energy increment of rubber sheet 
(1st load step). 








GLOBAL ENERGY INCREMENT ( INCH-LBS.) 


1st load step 



Figure 5.- Global energy increment of rubber sheet 
(1st, 2nd, 3rd load steps). 
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